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Lattice QCD allows us to study QCD phenomenology from first principles by using Monte 
Carlo techniques. Recent developments in both the computer technology and numerical algorithms 
have made possible lattice simulations with the correct number of fermion flavors in the vacuum 
polarization, which are essential for establishing direct connections between lattice simulations and 
the underlying low-energy QCD. However, the computational cost increases dramatically as one 
decreases the quark masses in the simulations towards the chiral limit. As such, current lattice 
simulations still work with quark masses heavier than their physical values, and extrapolations are 
necessary to obtain meaningful physical results from the simulations with heavy quark masses. 

Chiral perturbation theory (;^PT) is a low-energy effective theory which connects physical 
observables to quark masses in explicit functional forms, and is a useful tool to guide the extrap- 
olations for the lattice QCD simulations. Since it is based on the approximate chiral symmetry of 
QCD, it is important to have a chiral fermion formulation on the lattice in order to make direct use 
of the continuum ;t;PT for the sake of the extrapolations. The domain wall fermion (DWF) formu- 
lation is well-suited in this regard, since it preserves exact flavor symmetry, and chiral symmetry is 
only mildly broken. Its chiral symmetry breaking effect can be quantitatively described by a small 
additive mass shift called the residual mass, m^es- Recent work has shown ^] that, to do chiral 
extrapolations for domain wall fermions, the only modification to the continuum ;^PT is to replace 
the input quark mass by the sum of the input quark mass and wires, leaving the number of low energy 
constants unchanged, at least up to terms of ^(ma) which can be viewed as next-to-next-to-leading 
order (NNLO). This is in contrast to the cases of Wilson fermions or staggered fermions, where, 
at next-to-leading order, a few new low-energy constants need to be introduced to account for the 
chiral symmetry (Wilson) or flavor symmetry (staggered) breaking effects. 

One of the challenges of chiral extrapolations is that it still remains inconclusive what the 
radius of convergence is for ;^PT. Previous results of domain wall fermion simulations [|l|] have 
evidence that ;^PT at next-to-leadig order (NLO) is not sufficient to describe pion masses heavier 
than 400 MeV. One question to ask is, how light should the pion masses (or quark masses) be for 
;^PT to achieve the desired accuracy at NLO? In this work we present results for the pseudoscalar 
meson masses and decay constants from recent domain wall fermion simulations with 2-1-1 dynami- 
cal flavors on the 24^ x 64 lattices at a fixed lattice spacing of about 0. 1 fm. The partially quenched 
pion masses in these simulations are as light as 250 MeV, which gives us an opportunity to check if 
;^PT is consistent with the lattice data at this lighter mass range. The agreement between the lattice 
data and the predictions of in turn enables us to determine physical observables and the light 
quark masses with better controlled extrapolation errors than naive linear fits. In this proceedings 
we combine two talks given at the Lattice 2007 conference, and show our attempts to locate the 
mass range where ;tPT (SU(3)xSU(3) and SU(2)xSU(2)) is applicable, followed by the determi- 
nations of fj,;, /k and the physical light (up/down and strange) quark masses. For other physical 
results obtained from these configurations see ^ and references therein. 

1. Numerical Details 

The gauge configurations on the 24-' x 64 lattices were generated using the same parameters as 
the previous simulations on the 16^ x 32 lattices [Q]. Specifically, we used the Iwasaki gauge action 
with j8 = 2. 13. The extent of the fifth dimension was = 16, and the domain wall height was fixed 
to aMj = 1.8. The dynamical strange quark mass, arris = 0.04, was tuned to be approximately its 
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physical value, and four values of the light dynamical quark mass, ami, were used to allow for the 
extrapolations in the light quark mass limit. The rational hybrid Monte Carlo (RHMC) algorithm 
was applied to generate all the ensembles. The details of the implementation for the RHMC were 
reported in [^. The number of thermalized trajectories, in molecular dynamics time units, for the 
ami = 0.005,0.01,0.02 and 0.03 ensembles is 3600, 3600, 1760 and 1760, respectively. 

In order to make full use of the partially quenched ;^PT formulae [|6|], we calculated hadron 
correlators with the input valence quark masses am^ y E {0.001,0.005,0.01,0.02,0.03,0.04}. The 
lightest input quark mass turns out to be about 1/10 of the strange quark mass when the residual 
mass is properly included. Focus will be given to the two ensembles with lightest sea quark masses, 
ami = 0.005 and 0.01, as these smaller quark masses are more likely to be within the regime 
where NLO ;tPT has reasonable convergence. For these two ensembles, all the non-degenerate 
meson correlators were constructed from all the different combinations of the six valence quark 
masses using a Coulomb gauge fixed wall source (W) and either a wall or local sink (L) as part 
of our weak matrix element project [Q]. The quark propagators used in these measurements were 
obtained from the sum of quark propagators computed from periodic and anti-periodic boundary 
conditions, to eliminate the boundary effects from the backward-propagating states. For clarity 
we will denote these correlators as "W-P-i-A". Additionally, degenerate hadron correlators with a 
Coulomb gauge fixed 16^ box source and a 16^-box or local sink were also calculated on these 
ensembles. Note that correlators constructed from a box source and a box sink violate translational 
invariance, therefore zero-momentum projection can not be guaranteed. We thus summed over the 
correlators with all the possible choices for the box sink to achieve the zero-momentum projection. 
The correlators with a box source were found to have better overlap with the ground states of 
the baryons, and were used to extract the mass of the Q. baryon, which we will utilize to set 
the lattice scale for our simulations. The measurements were done on 90 gauge configurations 
on the ami = 0.005 and 0.01 ensembles, and about 45 on the am/ = 0.02 and 0.03 ensembles, 
with two adjacent measurements separated by 40 molecular dynamics time units. For each type 
of measurement we used two different source locations to reduce the fluctuations within the gauge 
configurations. The small number of measurements on each ensemble does not allow us to study 
the autocorrelation time reliably. However, we have checked that blocking the data in intervals of 
40 or 80 molecular dynamics time units does not change the statistical errors significantly, which 
is consistent with the study on the smaller volume [||]. Thus in the following analysis, we choose 
to block the data into intervals of 80 molecular dynamics time units for the ami = 0.005 and 0.01 
ensembles where 90 measurements are available, and into intervals of 40 molecular dynamics time 
units for the 0.02 and 0.03 ensembles, leaving approximately 45 jackknife samples for each of the 
four ensembles. 

2. Data Analysis 

2.1 The Residual Mass m^es and Axial Current Renormalization Za 

As the gauge coupling of this large volume simulation is identical to the 16^ x 32 simulations 
in [Q], we expect the residual mass mres and the axial current renormalization Za to be consistent 
with the results therein up to possible finite volume effects. The residual mass is determined from 
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the ratio [g] 

(E,J5"»')i"(0,0)> 

where J=,q is a point-split operator for domain wall fermions. Figure |T] shows the results oiR{t) at 
the four unitary points with anix = amy = ami. The horizontal lines represent the fit to a constant 
from f = 10 to 32 for each quark mass, determining am[^^{ami). The mass-independent residual 
mass is given by evaluating am[^^{ami) at amj = 0, and we have 

am,es = 0.00315(2). (2.2) 
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Figure 1: The ratio R{t) used in the determina- 
tion of the residual mass at the unitary points. 



Figure 2: Results for afp using different meth- 
ods on the ami = 0.005 ensemble. 



The axial current renormalization constant Za relates the local axial vector current 
q{x)x" y^y=,q{x) to the conserved axial vector current of domain wall fermions It is determined 
from the ^(a^) improved ratio as described in Ref. [^]. Similar to the residual mass, we compute 
the value of Za at each unitary quark mass, and extrapolate to the chiral limit at ami = —am^^s, 
obtaining 

Za = 0.7161(1). (2.3) 

The results for mres and Za are consistent with the previous results in the small volume, indicating 
no measurable finite volume effects are present for these quantities. 



2.2 Pseudoscalar Meson Masses and Decay Constants 

This section is devoted to the details of our fitting procedures to obtain the light pseudoscalar 
meson masses, mp. There are two types of interpolating operators which overlap with the pseu- 
doscalar meson state: P{x,y) = q{x)x''Y5q{y) and A{x,y) = g{x)z''Y4Y5q{y), where the quark fields 
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may have different smearings. To minimize systematic errors arising from different characteris- 
tics of the operators, we chose to fit all the available W-P+A correlators simultaneously to obtain 
a common mass and an amplitude for each correlator. Since the correlators are measured on the 
same gauge background, in principle we should take into account correlations among different 
correlators and different time slices of the same correlator. However, the fit range for a typical 
simultaneous fit is as large as 220 time slices. Having only 45 jackknife samples is not enough to 
resolve the covariance matrix. Thus the correlated fits cannot be carried out, and we restrict our- 
selves to uncorrelated fits. One caveat of the uncorrelated fits is that ;^^/d.o.f. from the uncorrelated 
fits does not follow the correct X'^ distribution, and do not reliably imply the goodness of the fits. 
The following five correlators were included in the simultaneous fits: 

(A^(OA^(0)), {P^{t)P^{0)), {A^{t)P^{0)),{P^{t)P^{0)), and {A^ {t)P^ {0)) , (2.4) 

where the superscripts indicate the smearing of the source or sink, with W being the Coulomb 
gauge fixed wall and L being the local operator. Each simultaneous fit gives a common mass amp, 
and one amplitude for each correlator, labelled as , £/p^ , £/^^ , and -s^^^ , respectively. 
There are five different ways to determine the pseudoscalar meson decay constant' , afp, using these 
amplitudes: 

^ ^''^ {amp)V .s/^p^ ^ {ampfV ^^^^ {amp)V i^^^ ^ 

^ {ampfV s^^p^ {amp)V ^ww2 ^^"^^ 

where V = {aL)^ is the spatial volume of the lattice. These ratios are calculated under a standard 
jackknife procedure to take into account correlations among different amplitudes. Note that not all 
of these methods are independent, but some of them may produce statistically more accurate results 
than the others due to different characteristics of the correlators. (II) and (IV) in fact come from 
the translation from pseudoscalar density to axial vector current using the axial Ward identity [^], 
hence the residual mass m^es is required. The results for afp of the ami = 0.005 ensemble from 
all of these different methods are shown in Fig. |2[ As we can see, they all give consistent results 
except that methods (II) and (IV) give slightly higher results than the rest at large masses, which 
may indicate different scaling errors resulting from the use of different correlators. In the following 
analysis, we use results from (III) since it gives the smallest statistical error. 

3. Chiral fits: SU(3) x SU(3) and SU(2) x SU(2) 

In this section we will discuss our attempts to fit the obtained meson masses and decay con- 
stants to formulae predicted by partially quenched chiral perturbation theory (PQ;^PT). (For similar 
fits for the kaon bag parameter Bk measured on the same lattice configurations see [l^.) Using 
PQ/^PT for three quark masses, corresponding to unquenched SU(3) x SU(3) ;tPT, up to NLO it 
will turn out that the data at our higher quark masses is not well described by the applicable formu- 



lae (Sect. p7l| ). Therefore, in Sect. 3.2 we will perform NLO SU(2) x SU(2) fits, dropping terms 
of order [mi/msY. 

'Our definition for the decay constant is such that the physical value of f-^ is about 130 MeV. 
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3.1 SU(3)xSU(3) Chiral Fits 

The most natural approach to fit our data from the Nf = 2 + \ ensembles is to use SU(3) x 
SU(3) xVT or its partially quenched variant, describing the dependence of the meson masses and 
decay constants on the two (in our case degenerate) light quark masses and the heavier strange 
quark mass by introducing chiral fit parameters to leading order (LO: Bq, fo) and next-to-leading 
order (NLO: L4 5 6,8), where the latter are commonly referred to as Gasser-Leutwyler parameters 
or low energy constants (LECs). From the general formulae given in the Nf = 2 + \ case has 
been worked out, see, e.g., |[l]]. 

When appying these fit forms to our data, we found that using (PQ);^PT to NLO does not 
describe our data well up to meson masses comparable to the kaon mass or — equivalently — up to 
an average quark mass of half the strange quark mass. We performed combined fits to a/^,, and 
{anixyY, meaning the decay constant or mass squared of a meson built from valence quarks with 
masses nix and niy, using the two ensembles with dynamical light quark masses of ami = 0.005 
and 0.01. A reasonable ;^^/d.o.f. could only be obtained by imposing a cut in the average valence 
quark mass of aniavg = {am^ + amy)/! < 0.01; fits with such a cut are shown in Fig. ^, while 
the fit parameters are given in Tab. |l|. There we conveniently quote the scale-dependent LECs 
at two commonly used chiral scales of = l.OGeV and 770 MeV. Also included in the table 
are phenomenological estimates for the LECs from JTl] ] and references therein. Our results show 
agreement with their NNLO fit values. In Fig. ^ we show fits with the cut chosen to be am^yg < 0.03. 
The fits miss almost all the data points inside the fitting range. Therefore, we conclude that NLO- 
XVY fits are not reliably applicable in a mass range up to the kaon mass. If one were to extract just 
the pion sector quantities, i.e., just the physical fj^, m^, and the physical light quark mass, from 
the fit results with the low mass cut, one still would include the terms proportional to the strange 
quark mass. Ideally, one would like to use ;^PT to guide the interpolation to the physical value of 
the latter. However, since we saw that at such a quark mass the fits deviate substantially from the 
data, this procedure has to be seen as an unsafe or at least questionable one. For the same reasons, 
a meaningful extraction of quantities in the kaon sector is impossible within this approach. 

One could try to extend the range of validity of ;^PT by going from NLO to NNLO. The 
complete formulae are available in the literature [12]. However, this would introduce much more 
LECs than the number of independent data points which are currently available to us. In addition, 
under these circumstances, we would not be able to establish whether such a NNLO fit was itself 
appropriate for this mass range. Were this kinematic region outside the domain of validity of ;tPT, 
such NNLO terms may not correctly describe our results. Instead we followed a different ansatz, 
namely to base the fit formulae just on the (approximate) chiral symmetry within the light quark 
doublet, as will be described in the next subsection. 



3.2 SU(2)xSU(2) Chiral Fits 

First, we will purely focus on the pion sector. By applying NLO SU(2) x SU(2) (PQ);^PT, 
where terms of order {mi/m^y have been dropped, the strange quark mass will not explicitly enter 
the fit formulae. The dynamical strange quark mass present in our simulations acts as a background 
field and is therefore implicitly contained in the SU(2) x SU(2) LECs. Of course, in that way we 
will not be able to correct for the difference between the dynamical strange quark mass value. 
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amy amy 




Figure 3: Combined SU(3) x SU(3) fits for the meson decay constants (left panels) and masses (right 
panels) at two different values for the light sea quark mass, valence mass cut amavg < 0.01. Points marked 
hy filled symbols were included in the fit, while those with open symbols were excluded. 



which was fixed during the generation of the gauge configuration (am., = 0.04), and its value at the 
physical point. As we shall see later on after extracting the physical amf^^^, this difference amounts 
to about 15 percent. 

The fit formulae for Nf = 2stre derived starting from too, as has been done, for instance, in 
[13]. Combined fits to af^y and {arrixy)^ from the two lightest ensembles with a mass cut of am^cvg < 
0.01 are shown in Fig. ^ whereas fitted parameters are included in Tab. |l[ Here we would like to 
point out, that — in contrast to SU(2) x SU(2) ;^PT — in the partially quenched theory the same set 
of LECs (instead of a reduced set due to operator redundancies) appears as in SU(3) x SU(3), since 
we have to distinguish between sea and valence quarks. This distinction gives rise to a different 
functional dependence of the considered quantities on the sea and valence quark masses. (More 
correctly speaking, we use a SU(4|2) x SU(4|2) chiral Lagrangian and not a SU(2) x SU(2) one.) 

Following this ansatz, as Fig. ^ suggests, we did not cure the problem of not being able to 
extend the fit range towards higher quark mass values. The important point is that our formulae 
do not contain any explicit dependence on the strange quark mass, whose physical value still lies 
outside the validity of the fit. The dependence on the background strange quark mass is implicitly 
contained in the LO and NLO fit parameters. One way (and in our opinion the most reliable one) 
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Table 1: Fitted parameters from different fits with a valence mass cut amavg < 0.01. For each fit the LECs 
are quoted at two different scales A;^. (Note: the value of Bq depends on the renormalization scheme like 
the quark masses: to obtain Bq, e.g., in the MS(2GeV) scheme, one has to divide the here quoted values by 



Z,„ (2GeV) from Sect. 4^.) Also included are the phenomenological estimates from [|ll|. Errors on Lg and 



L$ in Ql 1[] are added by quadrature to give the error on IL^ — L5. 
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Figure 4: Combined SU(3) x SU(3) fits for the meson decay constants (left panels) and masses (right 
panels) at two different values for the light sea quark mass, valence mass cut amavg < 0.03. Points marked 
by filled symbols were included in the fit, while those with open symbols were excluded. 



8 



Chiral Limit and Light Quark Masses in 2+1 Flavor Domain Wall QCD 



M.F. Lin and E.E. Scholz 



to estimate this dependence would be to repeat the same analysis on a second set of ensembles, 
generated at a different dynamical value for arris. 



To compare these fit results with the previous ones obtained in SU(3) x SU(3) ^PT (Sect. |3.1[ ), 
we use the formulae quoted in [|l^] (cf. also [H^]) to match the three flavor ;^PT to the two flavor 
case at LO. The results for aBo, afo and the low energy scales I3 4 (for a definition of the latter see 



[]14[]) are shown in Tab. g. From the fact that the converted SU(3) x SU(3) fit results almost agree 
with the SU(2) x SU(2) fit results, one may argue that the effect of a shghtiy too high strange quark 
mass may be neglected for quantities in the pion sector. 

Turning the attention now towards the kaon sector, we will have to incorporate the strange 
quark mass value. Since we already saw that NLO-PQ^PT fails to describe our data in the region 
of the kaon mass or even beyond, we decided to demand chiral symmetry properties only for 
the two light quarks. Analogously to the heavy-light chiral perturbation theory in the B-sector 
[ p^ p7| ] we propose to use SU(2) x SU(2) ;t;PT in the presence of K mesons with terms of order 
{mjt/niK)'^ being dropped at NLO. In other words, the K mesons are now not treated as pseudo- 
Goldstone bosons. Under such considerations, we give, in the following, the fit formulae for the 
decay constant and squared mass of a meson made from a light valence quark with a mass arrix 
and a heavier valence strange quark {anis). Here the dynamical light quark mass {ami) is taken into 
account as well, but the dynamical strange quark is viewed as a background field. (We followed 
the same ansatz to fit the kaon bag parameter Bk [0, p^.) 



1 



Xx + Xi, Xx + Xi , Xi - 2Zv , Xx 
log 0/ > N7 H log ■ 



(47ra/o)2L 2 2{aA^)^ 4 (aA;^)^ 

r,m. 111 "1 .,1 2 



(3.1) 

- aB™^<j 1 + f^xi + T^Zx \. (3.2) 

Here the mass parameter Xx = 2aBQ{amx + am^,.^) is used. The fit parameters a/o/f, aB^f^, ci 2, and 
d\^2 all carry a superscript rtis to indicate that these depend on the strange quark mass value. The 
parameters a/o and uBq are the same as the SU(2) x SU(2) ;^PT in the pure pion sector. Actually, 
in the following we fixed these to their values previously determined in the fits of the pure pion 
sector. 



We will use Eqs. (3.1, 3.2) to extrapolate the kaon decay constant and mass to the physical 
value of the light quark masses at a fixed value of the strange quark mass am^. Repeating this for 
different values of arrig allows us then to interpolate to the physical strange quark mass point as well. 
For the moment, since we only have data at one value for the dynamical strange quark mass, we 
can only vary the valence strange quark mass. For future runs, one should consider to have at least 
two sets of ensembles at different values of artis to allow for an interpolation between dynamical 
strange quark mass points. Finally, the fits were performed at arris = 0.03 and 0.04, using all the 
points with light (dynamical, valence) quark masses ami x < 0.01, i.e., the two ensembles with the 
lightest quark masses. Such fits at ants = 0.04 are shown in Fig. ^ where the diamonds indicate the 
extrapolations to the physical light quark mass, am^^^^, at the fixed anis. By interpolating between 
the results at these two values for the strange quark mass, we are able to extract afn and anin at a 
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Figure 5: Combined SU(2) x SU(2) fits for the meson decay constants {left panels) and masses {right 
panels) at two different values for the light sea quark mass, valence mass cut amavg < 0.01. Points marked 
hy filled symbols were included in the fit, while those with open symbols were excluded. 

Table 2: Comparison of converted SU(3) x SU(3) fit parameters with those from SU(2) x SU(2) fits. Low 
energy scales 1^ 4 are defined at A = 139MeV. 







aBo 


afo 


h 


h 


SU(3) 


X SU(3), conv. 


2.457(78) 


0.0661(18) 


2.87(28) 


4.10(05) 


SU(2) 


X SU(2) 


2.414(61) 


0.0665(21) 


3.13(33) 


4.43(14) 



determined physical strange quark mass, amf^'^^ , or, vice versa, use either the physical value (given 
the lattice spacing) of or niK to set ^^m^'^^^ 



4. Obtaining Physical Results 

First, we will discuss how the lattice scale was set and the points of physical quark masses 
were determined. In the remainder of this section we will utilize a non-perturbative renormalization 
scheme (RI/MOM) to obtain the quark masses in the MS scheme at /i = 2 GeV. 

4.1 Determination of amf^^^, amf^^^, c?"^ 

Given the reservations to use either the p -meson mass (width of the resonance) or the Sommer- 
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Figure 6: SU(2) x SU(2) fits for the kaon sector. Left panel: Kaon decay constant, right panel: Kaon mass 
for ams = 0.04. Points with filled symbols were included in the fit, while those with open symbols were 
excluded. 

scale (ten percent uncertainty due to phenomenological models) to set the lattice scale, here we will 
use the mass of the baryon, a state made out of three strange quarks, instead. One advantage of 
using this baryon mass, is that up to NLO in ;^PT it is free of logarithms containing the light quark 
masses [p^]. Therefore the extrapolation of the measured masses to the light physical mass can be 
easily performed using a linear ansatz without an uncertainty due to chiral logarithms. We used 
the measured masses [19] with anis = 0.03 and 0.04 extrapolated to the light physical masses 
(using the configurations with ami = 0.005 and 0.01) to interpolate to the value of the physical 
strange quark mass. 



The quark masses were obtained from the SU(2) x SU(2) fits described in Sect. ^2. For the 
light quark mass we solved for a pion mass of niji = 135.0MeV, corresponding to the physical 
uncharged pion mass, while for the strange quark mass the fit to the kaon mass was solved at 
mk = 495.7 MeV, which is the quadratically averaged neutral and charged kaon mass. 

Since these two determinations depend on each other (the lattice scale is needed to convert 
the input masses into lattice units, whereas the quark masses are needed for the extrapolation in 
the light and interpolation in the strange quark masses for the baryon mass), we performed these 
two steps iteratively, starting with an initial guess for the quark masses. After eight iterations no 
further relevant change in the parameters were observed. The final values for I /a, a, ami, amg can 
be found in Table |3| (including only the statistical error). 

Finally, with the knowledge of the values for the quark masses corresponding to their physical 
values, our chiral fits were used to extrapolate the meson decay constant to fji = 124. 1(3.6) MeV 
and interpolate to /k = 149.6(3.6) MeV (statistical error only). Compared to their experimentally 



observed values [gOp of 130.7(0. 1)(0.36) and 159.8(1. 4)(0.44) MeV, our values are about five or 
six percent too low, but our measured ratio {afK)/{afji) = 1.205(18) agrees within the uncertainty 
with the experimental value of 1.223(12), indicating possible scaling effects in our results. 

An interesting application of the latter result is to use it for the determination of the ratio 
|K.s of CKM-matrix elements, as has been pointed out in [pT|]. Using the input for the 

branching ratios Y{K — > }jLv{y)) and r(7r i^v{y)) plus radiative electroweak corrections from 
[20], we obtain |V„a = 0.2292(35) from our result for the decay constant ratio. This implies 
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Table 3: Determined lattice scale and spacing and unrenormalized quark masses {am^^^ = am^^'^ +flmres)- 

a-VGeV ^/fm am^^'' amf^' am^/'' amf^' 

1.729(28) 0.1141(18) -0.001847(58) 0.001300(58) 0.0343(16) 0.0375(16) 



\Vus\ = 0.2232(34), if = 0.97377(27) from super-allowed nuclear j3-decays ^ is taken into 
account. The quoted error combines both the errors from our determination of /k/ fn (statistical 
only) and the other input quantities. Here the main contribution comes from the decay constants, 
e.g., in the case of \Vus \ its contribution is 0.0033, whereas the other errors add up to 0.0005. 

4.2 Non-Pertubative Renormalization and Quark Masses 

The renormalization factor Z,„ = \/Zs needed to convert the extracted quark masses to the 
commonly used MS scheme at a scale of 2 GeV has been calculated (amongst others) using 
16^ X 32, Li = 16 DWF configurations with Nf = 2 + \ flavors [22]. (For details on the used 



configurations, cf. f^].) We first matched the bare lattice operators to the RI/MOM scheme using 



the non-perturbative Rome-Southampton technique []23|], followed by a perturbative matching to 
the MS scheme. Since DWF were used, we benefit from the controlled (small amount of) chiral 
symmetry breaking, resulting in &{a) improved operators/currents with reduced operator mixing. 

In particular, we calculated the renormalization factor in the regularization independent 
(Rl-)scheme according to 

Z^iap) = ^{ap)^iap)^, (4.1) 

where the first two factors were obtained from the renormalized amputated vertex functions A^*^" 
and A^^" (A^'^" = (Zv/Z^)Av), respectively, and the last factor was obtained by measuring the ap- 
propriate hadronic matrix element (see [Q]). The four loop matching from the RI to the renormal- 
ization group invariant (RGI-)scheme has been applied to extract Zm^^'^°'^(2GeV), which then 
was converted to the MS-scheme via three loop matching [^]. Finally, we get Z^^(2GeV) = 
1.656(48)(1 1), where the first error is the statistical one and the second one estimates the system- 
atics due to residual chiral symmetry breaking. The latter was obtained from the difference which 
arises if instead of A^ the combination (A^ + Ai/)/2 is used in the determination of Z„,. 
Using this result combined with the lattice spacing we obtain the quark masses via 



Z^^ (2 GeV) • ( 1 /a) • amf^' . (4.2) 



The physical light quark mass (which, in fact, is the average up- and down-quark mass) we measure 
is mi = 3.72(16)MeV, while for the strange quark mass we get a value of m^ = 107.3(4.5) MeV. 
(The quoted errors include the combined error from but only the statistical ones from other 
quantities.) This means we observe a quark mass ratio of m; : = 1 : 28.8(4). 

Conclusions & Outlook 

After realizing that fits to NLO ;tPT for three flavors are problematic up to the physical strange 
quark mass once sufficiently light quark masses have been reached, we found that using two flavor 
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XVT for the pionic sector is a much more rehable approach. It ehminates the questionable depen- 
dence of the pion mass and decay constant on the strange quark mass value explicitly. Nevertheless, 
by converting our three flavor ;^PT fit parameters to the two flavor case, a sufficient agreement be- 
tween the two approaches could be established. We quoted the fitted LO and NLO parameters for 
both the two and three flavor case. 

By only demanding chiral symmetry properties for the two hghtest quarks, we were able to 
apply SU(2) x SU(2) ;^PT to the kaon sector and successfully extracted the kaon mass and decay 
constant, despite the caveat that we had to include partially quenched strange quarks in that analysis 
because currently we are lacking data at a second value for the dynamical strange quark mass. 

By using the experimentally measured values for nij^, tuk, and m^j-, we were able to ex- 
tract the physical average light quark mass and strange quark mass, where for the conversion to 
the MS(2GeV) via the RI/MOM scheme a non-perturbative renormalization technique was used. 
The pion and kaon decay constants were extrapolated or interpolated to these quark mass values. 
We also derived the ratio of CKM-matrix elements |Vi„|/|V„d| from fx/ fn- For the moment, no 
estimates for systematic errors (except for Z,„) are given, which we shall do in a forthcoming pub- 



lication 1 25] 



Currently, we are running simulations at a larger lattice volume (32 x 64, = 16), where 



also a second value for the strange quark mass will be included. (For a status report see [26].) 
These simulations will allow us to estimate the size of finite volume errors as well as to interpolate 
between dynamical strange quark mass values, resulting in a more reliable result for the kaon 
sector. 
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